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Abstract 


The  method  of  determining  the  thickness  of  diffuse  phase 
boundary  with  density  profile  governed  by  equilibrium  condition?, 
is  proposed.  It  follows  the  well-known  procedure  of  analyzing  the 
deviations  from  Porod's  law.  Frrors  in  the  obtained  boundary  thick¬ 
ness,  due  to  the  statistical  scatter  in  the  intensity  data  and  to 
the  difficulty  of  separating  the  effect  of  density  fluctuations  with 
in  the  phases,  are  examined.  For  this  purpose,  scattering  curves 
are  synthesized  on  the  basis  of  a  well-defined  model  structure  with 
known  boundary  thicknesses.  These  synthesized  curves,  when  analyzed 
according  to  the  proposed  method,  yield  the  correct  boundary  thick¬ 
ness  under  favorable  conditions,  but  are  also  shown  to  lead  to  very 
erroneous  results  in  some  cases. 


I.  Introduction 


The  intensity  of  small-angle  X-ray  scattering  (SAXS)  fcv  iso¬ 
tropic  materials  having  Ideal  two  phase  structure  decreases  as  s  ’ 
at  large  s  Cs  =  2sinQ/l) ,  according  to  the  theory  by  Porod  (1951, 
1952a,  b) .  Deviations  from  this  Porod  law  can  occur  in  practice 
since  the  boundaries  between  the  phases  may  not  be  perfectly  sharp, 
and  the  electron  densities  within  the  phases  contain  local  fluctu¬ 
ations.  The  potential  utility  of  determining  the  phase  boundary 
thickness,  by  examination  of  the  deviation  of  observed  SAXS  inten¬ 
sity  from  Porod' s  law,  was  first  pointed  out  by  P.uland  (1971)  and 
was  since  examined  by  a  number  of  workers  (Vonk  1973;  Ruland  1974; 
Hashimoto,  Todo ,  Itoi  and  Kawai  1977;  Todo,  Hashimoto  and  Kawai 
1978;  Hashimoto,  Shibayama  and  Kawai  1980;  Koberstein,  Morra  and 
Stein  1980;  Roe,  Fishkis  and  Chang  1931). 

According  to  Ruland  (1971),  the  observed  intensity 
can  be  expressed  as 

IobsCs)  =  lidCs)H2(s)  +  lb  (1) 

where  I^(s)  is  the  intensity  due  to  ideal  two  phase  structure, 

H  (s)  is  the  factor  arising  from  diffuseness  of  the  boundaries, 
and  Ib  is  the  ’background’  intensity  due  to  local  electron-density 
fluctuation  within  the  phases.  If  the  electron-density  profile 
across  the  boundary  approximates  a  convolution  product  of  a  step 
function  and  a  gaussian  function,  H  (s)  takes  the  form 

K2  (s)  =  exp (-4t2c2s^ )  (gaussian)  (2) 

where  the  variance  a  of  the  gaussian  function  embodies  the  diffuse¬ 
ness  of  the  boundary.  Statistical-mechanical  studies  of  inhomo¬ 
geneous  liquid  systems  (Cafcn  and  Hilliard  1953)  ,  polymer  mixtures 
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tHelfand  and  Tagarai  1971;  Hopper  and  Uhlmann  1974;  Roe  1975)  and 
block  copolymers  (Helfand  1975)  showed  that  the  density  profile 
across  the  phase  boundary  attained  under  equilibrium  conditions 
is  given  by 

g(x)  =  tanh(irx//IT  cr)  (3) 

In  this  case  the  corresponding  H^(s)  function  is  (Ruland  1980;  Roe, 
Fishkis  and  Chang  1981) 

Cs)  =  (/ITn-as)  ^csch^  (/Terras)  (equilibrium)  (4) 

The  boundary  thickness  (or  the  value  of  a)  is  determined  by 
analysing  the  observed  intensity  in  the  Porod's  region  according 
to  eq.  (11,  and  for  this  purpose  a  suitable  graphical  or  numerical 
method  has  to  be  devised.  The  simplest  is  to  expand  H  (s)  in  a 
series  and  retain  the  first  two  terms  only.  Kober stein,  Morra  and 
Stein  (19801  showed  that  such  a  procedure  is  likely  to  lead  to  an 
underestimate  of  cr  when  the  analysis  is  performed  over  a  range  of 
s  in  which  truncation  of  higher  order  terms  is  not  justified.  They 
instead  proposed  a  procedure,  applicable  when  (s)  is  given  by 
ec.  (2) ,  which  is  valid  even  for  s  greatly  exceeding  1/c.  One  of 
the  purposes  of  the  present  work  is  to  present  a  similar  procedure 
which  is  applicable  wher,  H  (s)  is  given  by  eq.  (4)  instead  of  eq.  (2)  . 

The  main  practical  difficulty  in  utilizing  eq.  (1)  for  esti¬ 
mation  of  a  arises  from  uncertainties  introduced  by  the  errors  in 
the  observed  data.  The  effect  of  diffuse  boundaries,  represented 
by  H  (s) ,  introduces  a  negative  deviation  from  Porod's  law,  while 
the  'background'  1^  contributes  a  positive  deviation.  These  two 
effects  therefore  tend  to  cancel  each  other,  and  meaningful  sepa¬ 
ration  of  one  from  the  other  is  often  found  difficult  because  of 


the  statistical  error  in  the  data,  which  in  the  Porod's  region  is 
usually  appreciable.  It  is  desirable  to  have  an  estimate  of  the 
error  in  the  evaluated  o  value  arising  in  practice  as  a  result  of 
the  statistical  error  in  the  data .  Por  this  purpose  in  this  work 
a  series  of  'theoretical*  scattering  curves  based  on  a  model  struc¬ 
ture  are  constructed,  and  are  then  given  random  statistical  fluctu¬ 
ations  to  make  than  simulate  experimentally  determined  data.  These 
synthesized  scattering  curves  are  thereafter  analyzed  by  the  method 
proposed  for  the  determination  of  cr  values,  and  the  values  thus  ob¬ 
tained  are  compared  with  the  values  initially  chosen  for  the  syn¬ 
thesized  curves. 


II.  Procedure  Applicable  to  Equilibrium  Boundary  Profile 

When  SAXS  data  are  collected  with,  instruments  having  slit 
collimation,  the  effect  of  the  slit-smearing  has  to  be  properly 
taken  into  account.  Because  of  the  possible  amplification  of  the 
error  in  the  data  in  the  process  of  desmearing,  the  procedure  for 
the  determination  of  the  boundary  thickness  should  preferably 
make  use  of  the  smeared  data  directly.  The  slit-smeared  inten¬ 
sity  ICs)  is  given  by 

ItS)  =  2J^  VI ^  tu)  Iobs C/si+u2)  du  (5) 

where  W^tu)  is  the  slit-length  weighting  function.  In  what  follows, 
we  will  consider  only  the  cases  in  which  the  inr inite-slit  approxi¬ 
mation  is  valid.  W^ tu)  is  then  a  constant  independent  of  u. 

The  Intensity  of  Ideal  two  phase  structure  in  the  Porod's 
region  is  given  by 

Iidts)  =  (Ap)2S/8tt3s4  (6) 

where  Ap  is  the  difference  in  the  electron  density  between  the  two 
phases  and  S  is  the  total  phase  boundary  area.  For  the  density  pro¬ 
file  across  the  phase  boundary  determined  by  equilibrium  conditions, 
K2 ts)  is  given  by  eq.  (4),  and  therefore  the  slit-smeared  intensity 
in  the  Porod's  region  can  be  written  as 

ICs)  =  Agtx)  +  Ib  C7) 

where  Ib  is  the  smeared  'background'  intensity,  and 

A  =  36irSCAp)2W^(0)a4  C8) 

x  =  /12-rras  (9) 
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2  "> 
g  Cxi  =  csch  (x)/x^ 

(101 

gCx)  =  P  g (/x2j-u2)  du 

(11) 

The  integration  indicated  in  ec.  (11)  could  not  be  evaluated 
analytically,  and  in  this  work  g(x)  was  calculated  numerically . 

For  x  between  0.8.  and  6.0,  g(x)  is  approximated  well  by 

gCx)  *  3.3x“1,8e~1,9*  (12) 

The  root-mean-square  error  of  this  approximation  over  the  above- 

mentioned  range  of  x  is  2%  and  the  maximum  error  at  the  two  ends  of 

1  8 

the  range  is  4%.  Fig.  1  shows  the  plot  of  ln[x  *  g(x)J  against  x, 

illustrating  the  good  linearity  in  the  indicated  range  of  x. 

Evaluation  of  a  from  the  observed  intensity  I(s)  would  then 

^  18 

proceed  as  follows.  A  plot  of  ln[{I(s)  -  I^ls  ’  ]  against  s  is 
constructed,  and  the  slope  of  the  linear  portion  of  the  plot  is 
determined.  The  value  of  a  is  then  given  by 

a  -  -slope/  (1. 9/Hi)  (12) 

Alternatively,  one  could  evaluate  the  effective  thickness  t  by 

t  ■  -slope  (2/1 . 9tt2)  (14) 

where  t,  defined  as  the  reciprocal  of  the  density  gradient  at  the 
center  of  the  boundary  layer,  is  related  to  c  by 

t  ■  (4/T /t)c.  (15) 


After  a  has  been  evaluated,  one  may  calculate  x  by  eq.  (9) 
to  see  whether  the  approximation  (12)  is  justified  in  the  range 
of  s  at  which  the  slope  was  evaluated.  In  the  (unlikely)  event 


that  the  x  range  starts  at  a  value  smaller  than  0.8,  the  following 
alternative  approximation  may  be  employed. 


~  ,  .  _  _  *  -2.5  — 1 . 4x 

gCx)  -  2, Ox  e  Cl 6) 

Eq.  Cl 6)  is  valid  for  x  between  0.3  and  2.1,  with  the  root-mean- 
•square  error  equal  to  2%  and  the  maximum  error  equal  to  4%.  The 
required  plot  is  then  In [{Its}  -  Ifc}s  *  ]  against  s. 

III.  Synthesized  Scattering  Curves 

It  would  be  a  difficult  task  to  derive  an  analytical  expression, 
in  general  terms,  for  the  propagation  of  error  from  the  initial  data 
to  the  final  a  value  obtained.  Instead,  we  take  a  simpler  approach 
of  synthesizing  a  model  scattering  curve  and  then  applying  the  pro¬ 
posed  procedure  to  it  to  see  how  closely  the  known  c  value  is  re¬ 
covered.  In  constructing  the  model  scattering  curves,  effort  is 
given  to  make  them  as  realistic  as  possible.  The  physical  constants 
required  for  their  calculation  are  given  values  typically  found  in 
real  systems  of  interest. 

Two  types  of  scattering  curves  are  constructed.  In  the  first, 
the  model  structure  assumed  is  one  in  which  spherical  particles  of 
one  phase  are  imbedded  in  a  continuous  matrix  of  another  phase,  with 
the  boundary  between  the  phases  given  a  predefined  diffuseness.  The 
second  corresponds  to  a  scattering  from  a  single-phase  two-components 
system.  This  is  to  examine  what  erroneous  a  values  are  likely  to  be 
obtained  when  the  system  actually  contains  no  two  phase  structure 
at  all. 

A.  Structure  containing  spherical  particles 

Construction  of  the  scattering  curves  proceeds  in  three  stages: 
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deriving  the  scattering  curve  for  pin-hole  collimation,  applying 
the  slit  smearing  and  adding  random  statistical  error. 

For  an  ideal  two  phase  system  containing  spherical  particles, 
of  radius  a  distributed  according  to  a  distribution  function  w(a), 
the  scattered  intensity  per  unit  volume  of  the  sample  is 

Iid(s)  =  (Ap)2  f1J(j>2(s,a)w(a)da/J’(4/3)^a3wCa)da  (17) 

2 

where  f^  is  the  volume  fraction  of  the  spheres,  and  $  (s,a)  repre¬ 
sents  the  scattering  from  a  sphere  of  radius  a,  given  (Guinier  1963) 
by 

4>ts,a)  =  C4/3)  TTa3u  C2Tisa)  CIS ) 

Utx)  =  (3/x^)  (sinx  -  xcosx)  (19) 

For  w(a) ,  we  take  a  triangular  distribution  rising  linearly  from 
aQ-Aa  to  aD  and  then  falling  linearly  from  aQ  to  aQ+Aa.  (This  is 
used  here  in  preference  to  the  usual  gaussian  because  it  imposes 
an  upper  and  lower  limit  to  a.)  Eq.  (17)  represents  the  sum  of 
independent  scattering  by  the  spheres  and  ignores  the  effect  of 
inter-particle  interference,  but  this  is  justified  because  in  the 
Porod's  region  only  the  intra-particle  scattering  makes  recog¬ 
nizable  contribution. 

The  pin-hole  collimation  scattering  curve  Its)  is  then  ob- 
tained,  according  to  eq.  (1),  by  multiplying  Ij^(s)  with  H  (s)  to 
account  for  the  diffuse  boundaries  and  by  adding  for  the  'back¬ 
ground.  *  Eq.  (4),  valid  for  the  boundary  profile  obtainable  under 
equilibr ium  conditions,  is  used  for  H  ts) .  1^  depends  weakly  on  s, 

and  its  value  extrapolated  to  s=0,  for  a  pure  amorphous  material, 
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is  given  CGuinier  1963)  by 

IbCO)  =  p2k T  scT  C201 

where  p  is  the  electron  density  and  <T  the  isothermal  compress¬ 
ibility. 

The  following  numerical  values  are  used  in  the  computation: 
ip  =  4  x  10"2  electrons/A^ ,  =  0.25,  aQ  =  1Q0A,  ia  =  20A,  I^CO)  = 

0.3  electrons^/A^ .  These  are  the  values  which  would  be  found  for 
example  with  styrene-butadiene  block  copolymers.  Three  different 
curves  with  boundary  thickness  t  equal  to  5,  20  and  40A  are  calcu¬ 
lated. 

The  slit  smearing  effect  is  introduced,  by  means  of  eq.  (5) , 
into  the  pin-hole  collimation  intensity  I(s)  thus  calculated.  When 
the  infinite-slit  approximation  is  taken,  W^Cu)  is  a  constant,  which 
is  here  set  arbitrarily  equal  to  0.5A.  For  the  'background1  I.  ,  the 
infinite-slit  approximation  cannot  be  used.  1^  and  1^  show  about 
the  same  weak  dependence  on  s,  but  their  relative  magnitudes  depend 
on  the  actual  slit  weighting  function  (u) .  For  the  slit  geometry 
typically  encountered  in  Kratky  cameras,  I j./ 1 j,  *  0.05  (when  W»  (0)  is 
equated  to  0.5A).  Fmpirically,  the  angular  dependence  of  was 
found  (Vonk  1973,  Todo,  Hashimoto  and  Kawai  1973,  Hashimoto, 
Shibayama  and  Kawai  1980)  to  obey  the  following  relation  fairly  well 

Ifc  =  a  fcs11  (21) 

with  n  an  even  integer.  The  following  values  of  the  constants  are 

used  in  this  work  for  the  numerical  calculation:  a  =  0.015 

2  ,  2 
electron  /A",  b  =  600  electron  A  and  n  =  4. 

Finally,  to  simulate  the  counting  statistical  error,  the  exact 


-3- 


Its)  values  calculated  at  discrete  s  values  are  replaced  by  random 
numbers  generated  in  accordance  with  gaussian  statistics  of  the 
appropriate  mean  and  variance.  In  the  preset- time  mode  of  data 
collection,  or  in  the  use  of  position-sensitive  detectors,  the 
variance  of  the  error  statistics  is  proportional  to  the  intensity 
at  each  data  point.  The  data  points  in  the  three  constructed  curves 
are  given  statistical  errors  accordingly,  with  their  magnitudes  ad¬ 
justed  to  make  the  standard  deviation  at  the  lowest  intensity  point 
equal  to  1%.  The  resulting  three  curves,  for  t  =  5,  20  and  40A,  are 
designated  SI,  S2  and  S3  for  later  reference,  and  are  shown  in  Fig. 

2. 

B,  Single-phase  structure 

Single  phase  solutions,  containing  polymer  molecules,  can  ex¬ 
hibit  appreciable  electron-density  inhomogeneity  due  to  concentra¬ 
tion  fluctuation  and  give  rise  to  scattering  of  X-ray.  The  intensity 
can  be  especially  pronounced  when  the  structure  is  on  the  verge  of 
phase  transition.  As  a  model  scattering  curve  from  such  a  system, 
we  take  the  theoretical  expression  derived  by  Leibler  (1980)  for  a 
homogeneous  block  copolymer  system  above  the  transition  temperature 
of  domain  formation.  The  intensity  of  scattering  given  by  eq.  (IV-5) 
of  his  work  depends  on  the  parameters  x  (the  interaction  parameter) , 

R  Cthe  radius  of  gyration) ,  N  (the  number  of  segments  per  molecule) 
and  f  Cthe  fraction  of  component  1) .  To  express  the  intensity  in 
electron  unit  per  unit  volume,  Leibler* s  expression  has  to  be  multi¬ 
plied  by  CAp)^Vu  where  Ap  is  the  electron  density  difference  between 
the  components  and  Vu  the  volume  of  a  polymer  segment.  The  follow¬ 
ing  numerical  values  for  these  parameters  are  chosen  for  computation: 
f  =  0.25,  R  =  100A,  N  =  10^,  Vu  =  10  A*'  and  Ap  =  0.04  electrons/^. 
Two  curves  are  constructed,  one  for  XN  =  17  and  the  other  for  x  =  0. 


The  second  corresponds  to  a  very  high  temperature  and  the  first  to 
a  temperature  just  above  the  transition  temperature. 

Slit-smearing  integration  is  then  performed  on  these  two  curves, 
the  'background'  is  added,  and  the  statistical  errors  are  intro¬ 
duced,  all  in  the  same  way  as  described  in  section  A  above.  Fig.  3 

shows  the  two  curves  finally  obtained.  They  will  be  referred  to  as 
* 

LI tXN  =  17)  and  L2 U  ■  0) . 

XV.  Error  in  the  Evaluated  Boundary  Thickness 

The  method  of  determining  the  boundary  thickness,  described  in 
Section  II,  is  applied  here  to  the  synthesized  scattering  curves  to 
see  how  well  the  original  thickness  values  are  recovered. 

A.  Effect  of  statistical  error 

~  -  1*8 

Fig.  4  shows  the  plot  of  In  [Cl  -  Ij>)  s  *  ]  vs.  s  for  the  three 
synthesized  curves  Si,  S2  and  S3.  The  'background'  Ij,  subtracted 
is  the  exact  one  which  was  initially  added  in  the  synthesis,  i.e., 
that  given  by  eq.  (21)  with  a  *  0.015,  b  =  600  and  n  =  4.  The  addi¬ 
tional  error  which  arises  when  a  slightly  erroneous  Iv,  is  subtracted 
is  to  be  examined  later.  Here  the  uncertainty  in  the  value  of  t 
Cor  a)  determined  arises  solely  because  of  the  large  scatter  in  the 
individual  data  points  which  becomes  aggravated  after  the  exact  1^ 
is  subtracted.  The  angular  range  of  data  for  which  the  slope  can 
be  evaluated  is  very  limited,  especially  in  the  case  of.  S3.  The 
data  in  the  low  angle  region  deviates  from  Porod's  law  because  of 
the  interference  effect  reflecting  the  shape  of  the  particles,  while 
the  high  angle  region  cannot  be  used  because  of  the  excessive  scatter 
of  the  points.  The  arrows  in  Fig.  4  indicate  the  range  of  s  used 
for  the  evaluation  of  the  slope.  The  boundary  thickness  calculated 
from  the  slope  by  ec.  C14)  is  3 ,  19  and  30  A  for  SI,  S2  and  S3 
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respectively.  These  values  are  considerably  different  from  the 
expected  values  5,  20,  and  40A.  For  curve  SI,  the  angular  range 
indicated  in  Fig.  1  corresponds  to  X  =  0.7 -1.3  when  s  is  converted 
to  x  by  eq.  (.9)  with  aC=  irt/A/J)  =  3.6.  The  lower  end  of  the  x 
range  thus  extends  slightly  beyond  the  validity  of  the  approximation 
given  by  eq.  (.12)  ,  and  this  partly  explains  why  the  t  value  obtained 
is  too  large.  For  S2  and  S3  the  error  in  t  probably  arises  simply 
from  the  large  uncertainty  in  the  determination  of  the  slope. 

To  see  whether  a  substantially  reduced  statistical  error  in 
the  data  can  lead  to  a  more  accurate  value  of  the  boundary  thick¬ 
ness,  additional  scattering  curves  similar  to  SI,  S2  and  S3  are 
synthesized  with  0.1%,  instead  of  1%,  random  counting  error.  The 
plots  according  to  eq.  (12)  for  these  new  curves  are  given  in  Fig. 

5.  Here  the  ranges  of  s  at  which  the  slope  can  be  evaluated  are 
wider  and  extend  to  larger  angles.  The  evaluated  t  values  are  6, 

19  and  36A,  which  approach  the  true  values  much  more  closely.  Even 
here,  however,  it  would  be  difficult  to  claim  an  accuracy  better 
than  10%  in  the  determination  of  t. 

B.  Effect  of  erroneous  background  subtraction 

In  constructing  the  plots  shown  in  Figs.  4  and  5,  the  same 
'background'  as  originally  added  in  the  synthesis  of  the  scattering 
curves  was  subtracted.  In  practice,  however,  the  'background'  has 
to  be  estimated  from  analysis  of  high  angle  data.  Error  in  the 
estimated  background  can  arise  not  only  because  of  the  scatter  in 
the  data,  but  also  because  the  valid  analytical  form  for  the  back¬ 
ground  intensity  curve  is  usually  not  known.  Several  workers  (Vonk 
1973;  Todo,  Hashiraoto  and  Kawai  1978;  Hashiraoto,  Shibayama  and 
Kawai  1980;  Roe,  Fishkis  and  Chang  1981)  found  the  empirical  eq.  (21) 
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to  fit  the  observed  background  fairly  well,  while  P.athje  and  Ruland 
(197  61  and  Wiegand  and  Ruland  (1979)  used  the  exponential  form 

rb  =  a  expCbs^}.  (.22) 

Among  the  four  possible  analytical  forms,  i.e.,  eq.  (.21)  with  n  =  2, 

4  and  6  and  ec.  (22) , 'more  than  one  may  give  apparently  satisfactory 
fit  to  the  data  within  the  error.  The  somewhat  arbitrary  choice 
of  the  analytical  form  can  then  influence  the  value  of  the  boundary 
thickness  finally  evaluated. 

In  the  synthesis  of  curves  SI,  S2  and  S3,  the  background  was 
added  by  eq.  (21)  with  n  =  4.  Fig.  6  gives  the  plot  of  these  scat¬ 
tering  curves  against  s6  as  the  abscissa,  and  shows  that  eq.  (21) 
with  n  =  6  is  a  fairly  satisfactory  approximation.  Here  the  data 
extends  up  to  s  s  0.05.  After  the  backgrounds  evaluated  according 

to  the  straight  lines  drawn  in  Fig.  6  are  subtracted  from  SI,  S2 

-  ^  18 

and  S3,  the  plots  of  ln[ (I  -  Ib) s  *  ]  against  s  are  constructed  in 
Fig.  7.  The  boundary  thicknesses  evaluated  from  the  slopes  are  12, 

22  and  49A,  respectively  all  substantially  exceeding  the  true  values 
5,  20  and  40A. 

In  Fig.  8  curves  SI,  S2  and  S3  are  replotted  to  show  that  eq. 

(22)  is  also  a  reasonably  good  approximation.  Here  again  s  extends 
up  to  0.05.  After  subtraction  of  the  background  given  by  the  straight 
lines  in  Fig.  8,  the  plots  to  evaluate  the  boundary  thickness  are 
constructed  in  Fig.  9,  from  which  extremely  erroneous  t  values  of  8, 

11  and  9A  are  obtained  for  SI,  S2  and  S3,  respectively. 

There  are  some  latitudes  in  the  values  of  a  and  b  which  are 
determined  from  such  plots  as  given  in  Figs.  6  and  3.  The  bound¬ 
ary  thickness  values  finally  evaluated  are  less  sensitive  to  the 
uncertainties  in  these  constants,  and  are  found  to  vary  less  than 
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10%  in  most  cases.  This  fact  probahly  led  to  the  occasional  state¬ 
ments  in  the  literature  that  the  boundary  thickness  values  deter¬ 
mined  were  not  very  sensitive  to  the  method  of  background  subtrac¬ 
tion. 

C.  Illusory  boundary  thickness  in  a  single  phase  structure 

An  experimental  scattering  curve  from  a  specimen  which  does 
not  contain  a  two  phase  structure  may  nevertheless  sometimes  yield 
an  illusory  boundary  thickness  if  the  procedure  is  blindly  applied 
to  it.  The  scattering  curves  synthesized  in  Section  IIIB,  on  the 
basis  of  the  theory  by  Leibler  (1980)  of  a  homogeneous  (single-phase) 
block  copolymer,  are  utilized  here  to  examine  this  possibility. 

Pig.  10  shows  the  two  curves  Ll  and  L2  plotted  against  s4,  to 
evaluate  the  background  in  accordance  with  ec.  (21)  with  n  =  4.  The 
values  of  the  constants  a  and  b  in  eq.  (21)  obtained  from  Fig.  10 
are  0.0196  and  280,  respectively,  which  are  appreciably  different 
from  0.015  and  600  used  in  the  synthesis  of  these  curves  initially. 
The  component  of  scattering  due  to  concentration  fluctuation  decays 
rather  slowly  with  s,  and  some  of  it  is  erroneously  included  in 
the  background  evaluated  in  the  above  manner .  The  erroneous  back¬ 
ground  is  subtracted  from  Ll  and  L2  and  Fig.  11  shows  the  plot  con¬ 
structed  for  the  determination  of  the  boundary  thickness  in  the 
usual  manner.  Quite  possibly  one  could  be  misled  to  taking  the 
slope  as  given  by  the  straight  line  there  and  come  to  the  conclusion 
that  t  is  equal  to  around  10A. 

V.  Discussion 

The  numerical  examples  given  in  the  previous  section  demon¬ 
strate  that  under  favorable  conditions  the  boundary  thickness  t  can 
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be  evaluated  correctly.  They  also  point  out  the  need  to  guard  against 
the  possibility  of  obtaining  seriously  erroneous  results. 


As  discussed  in  the  introduction,  the  basic  problem  is  that  the 
two  types  of  deviations  from  Porod's  law  due  to  the  diffuse  boundary 
and  the  density  fluctuation  produce  opposite  effects,  which  can  be 
separated  from  each  other  only  with  difficulty.  As  the  numerical 
examples  in  Section  IIIB  illustrate,  assumption  of  wrong  analytical 
forms  for  the  background  intensity  can  lead  to  grossly  misleading 
values  of  the  boundary  thickness.  The  danger  of  the  wrong  back¬ 
ground  equation  can  be  alleviated  by  extending  the  collection  of  the 
data  to  higher  angles.  This  may  not  eliminate  it  completely,  how¬ 
ever,  since  there  is  never  a  guarantee  that  an  empirical  equation 
giving  a  good  fit  at  high  angles  is  indeed  a  valid  one  in  the  low 
angle  region  of  importance  for  the  boundary  thickness  determination. 
If  each  of  the  components  forming  the  two  phase  structure  can  be 
obtained  in  pure  form,  more  reliable  background  can  be  obtained  by 
measuring  the  scattering  from  the  pure  materials  separately  and  cal¬ 
culating  the  weighted  average  of  the  two. 

The  statistical  error  of  counting  present  in  the  data  makes 
the  evaluation  of  the  thickness  doubly  difficult.  First,  it  makes 
the  separation  of  the  background  uncertain.  Second,  even  when  the 
background  can  be  subtracted  correctly,  the  scatter  in  the  data 
reduces  the  precision  of  the  obtained  t  value  markedly.  The  ex¬ 
amples  given  in  Section  IIIA  show  that  a  statistical  error  cor¬ 
responding  to  1%  standard  deviation  in  the  background  region  is 
barely  tolerable,  if  an  accuracy  in  t  better  than  20%  is  desired. 

One  should  probably  aim  at  achieving  a  statistical  error  of  less 
than  0.1%  in  the  background  —  an  aim  which  may  not  easily  be 


c 

attained  in  practice.  To  accumulate  10  counts  per  channel  in  the 
background  with  amorphous  polymers,  it  may  take  several  days  with 
a  Kratky  camera  fitted  with  a  one-dimensional  position  sensitive 
detector  and  operating  on  a  rotating  anode  tube.  Accumulation  of 
such  high  counts,  for  only  a  few  samples,  may  also  alter  the  detec¬ 
tion  efficiency  characteristics  of  the  position-sensitive  detector 
(Hoe,  Chang,  Fishkis  and  Curro  19fll)  noticably.  An  important  factor 
affecting  the  error  in  t  is  the  relative  magnitudes  of  the  Porod  Law 
intensity  and  the  background  intensity.  Collimation  with  a  long  slit 
will  increase  the  latter  more  than  the  former,  and  therefore  short 
slit  or  pin-hole  collimation  would  be  advantageous  under  otherwise 
comparable  conditions . 

The  conclusions  derived  in  this  work  are  based  on  numerical  re¬ 
sults  obtained  with  synthesized  scattering  curves.  Although  they 
are  constructed  to  model  typical  experimental  curves  in  realistic 
detail,  there  would  of  course  be  experimental  results  which  differ 
greatly  from  them.  In  such  cases  the  assessment  of  the  error  given 
in  this  work  may  not  be  exactly  valid.  One  may  then  be  required  to 
undertake  an  assessment  of  the  error  in  the  evaluated  boundary  thick¬ 
ness  by  repeating  what  is  described  in  this  work,  that  is,  by  syn¬ 
thesizing  a  scattering  curve  embodying  the  structure  under  study  and 
the  slit  geometry  used  and  by  examining  how  well  the  boundary  thick¬ 
ness  built  into  the  scattering  curve  can  be  recovered. 
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Legend  to  Figures 


Fig.  1.  The  plot  illustrates  that  eq.  C12]  approximates  g(x)  well 
for  x  between  0.8  and  6.0.  Here  g(x),  given  by  eqs.  CIO) 
and  0.1} /  is  the  function  representing  the  Porod's  law  in¬ 
tensity  for  structures  having  diffuse  boundary  in  the  case 
of  infinite-slit  geometry. 

Fig.  2.  The  plots  show  the  synthesized  scattering  curves  which  simu¬ 
late  the  SAXS  patterns  obtainable  with  infinite-slit  colli- 
mation  from  samples  containing  spherical  domains  (with  a 
triangular  distribution  of  the  radius  between  80  and  12 OA) 
imbedded  in  a  continuous  matrix  of  another  phase.  The  back- 

4 

ground  intensity  in  the  form  a+bs  is  assumed.  Curves  SI, 
S2,  and  S3  are  for  the  boundary  thicknesses  of  5,  20  and 
40A  respectively.  The  statistical  error  added  corresponds 
to  1%  standard  deviation  at  the  point  of  lowest  intensity. 

Fig.  3.  The  plots  show  the  synthesized  scattering  curves,  which 
simulate  the  SAXS  pattern  obtainable  with  infinite  slit 
collimation  from  a  sample  containing  a  homogeneous,  single¬ 
phase  block  copolymer.  The  theoretical  ecuation  derived  by 

4 

Leibler  (1980)  is  utilized,  but  the  background  (a+bs  ) , 
slit-smearing  and  statistical  errors  are  also  added.  Ll 
is  for  XN  =  17  and  L2  for  XN  =  0 . 


Fig.  4.  Curves  SI,  S2  and  S3  are  plotted  in  accordance  with  the 

method  described  in  Section  II  for  evaluation  of  boundary 
thickness.  (Por  the  sake  of  legibility,  the  plots  for  S2 
and  S3  are  displaced  downward  by  multiplication  of  0.1  and 
0.01  respectively.)  The  eucact  background,  as  was  initially 
added  in  the  construction  of  the  synthesized  ctn:es ,  was 


subtracted  from  them  before  making  the  plot.  The  range  of 
s  at  which,  the  slope  is  evaluated  is  indicated  5v  arrows. 

The  thickness  t  obtained  is  8 ,  19  and  3  0A  for  Si,  S2  and 
S3,  respectively. 

Fig.  5.  These  plots  are  similar  to  those  shown  in  Fig.  4,  except 

that  the  curves  Si”’,  32^,  and  S3"”  shown  here  are  ones  con¬ 
structed  with  statistical  errors  corresponding  to  0.1%  at 
the  lowest  intensity  data  point.  The  much  smaller  scatter 
enables  the  slope  to  be  evaluated  in  wider  ranges  of  s  ex¬ 
tending  into  larger  angles.  The  thickness  values  evaluated 
6,  19  and  36A  now  agree  better  with  the  true  values. 

Fig.  6.  The  plots  illustrate  that  a  reasonably  good  approximation  to 
the  background  can  be  obtained  even  when  a  wrong  analytical 
equation  ta+bs^,  instead  of  the  correct  one  a+bs^)  is  assumed. 
The  plots  for  S2  and  S3  are  displaced  downward  for  legibility 
by  subtraction  of  a  constant  term  0.01  and  0.02  respectively. 

Fig.  7.  The  erroneous  backgrounds,  evaluated  according  to  Fig.  6, 
were  subtracted  from  the  scattering  curves  SI,  S 2  and  S3, 
before  constructing  the  plot  to  evaluate  the  boundary  thick¬ 
ness.  (.The  plots  for  S2  and  S3  are  displaced  downward  for 
legibility.)  The  thickness  t,  evaluated  from  the  slopes 
indicated  by  the  straight  lines  drawn,  are  12,  22  and  49A 
for  the  three  curves  respectively,  all  larger  than  the 
true  values. 

Fig.  8.  The  plots  illustrate  that  the  background,  originally  synthe- 

4 

sized  according  to  a+bs  ,  can  be  fitted  fairly  well  by 
2 

a  exp  (bs  )  .  (The  plots  for  S2  and  S3  are  displaced  dcwnv'ard 
for  legibility  by  multiplying  by  0.5  and  0.25  respectively.) 
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Pig.  9.  The  erroneous  backgrounds ,  evaluated  according  to  Pig.  8, 


were  subtracted  from  the  scattering  curves  SI,  S2  and  S3, 
before  constructing  the  plot  to  evaluate  the  boundary 
thickness.  CThe  plots  for  S2  and  S3  are  displaced  down¬ 
ward  for  legibility! .  The  thickness  t,  evaluated  from  the 
slopes  indicated  by  the  straight  lines  drawn,  are  8,  11 
and  9A  respectively,  the  latter* two  being  grossly  in  error. 

Pig.  10.  The  single-phase  scattering  curves  Ll  and  L2  were  originally 
synthesized  with  the  background  0.015+600s4.  The  plots 
shown  here  illustrate  that,  despite  the  correct  equation 
with  a  4th  order  term  assumed  for  the  background,  the  best 
fit  obtained  differs  greatly  and  is  given  by  0.020+280s4. 
CThe  plot  for  L2  is  displaced  downward  by  0.01  for  legi¬ 
bility.  ) 

Fig.  11.  The  single-phase  scattering  curves  Ll  and  L2,  from  which 
the  erroneous  backgrounds  evaluated  in  Fig.  10  were  sub¬ 
tracted,  are  plotted  according  to  the  method  to  be  used 
for  evaluation  of  the  boundary  thickness.  CThe  plot  for 
L2  is  displaced  downward  by  multiplication  of  a  factor  0.1.) 
The  downward  trend  of  the  data  points  at  high  angles,  pro¬ 
duced  by  subtraction  of  the  incorrect  background,  can  erro¬ 
neously  be  interpreted  to  mean  the  presence  of  phase  bound¬ 
ary  of  thickness  11A. 


-20- 


(l-lb)s 


|S  (  q'|  —  | ) 


Fig.5 


Fig.  8 


(l-lb)s' 


r 


Fig.  9 


SP472-3/A1 


472  :  GAN* :  716  :  dec 
78u4?2-6C$ 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  GEN 


No . 

Copies 

Office  of  Naval  Research 

Attn:  Code  472 

300  North  Quincy  Street 

Arlington,  Virginia  22217  2 

0NR  Branch  Office 
Attn:  Dr.  George  Sandoz 
536  S.  Clark  Street 

Chicago,  Illinois  60605  1 

Office 

Attn:  S^Theof^fcDept . 

^<^YorkThTew  Yor^lTStro^  1 

ONR  Western  Regional  Office 
1030  East  Green  Street 

Pasadena,  California  91106  1 

ONR  Eastern/Central  Regional  Office 
Attn:  Dr.  1.  II.  Peebles 

Building  114,  Section  D 
666  Summer  Street 

Boston,  Massachusetts  02210  1 

Director,  Naval  Research  Laboratory 
Attn:  Code  6100 

Washington,  D.C.  20390  1 

The  Assistant  Secretary 
of  the  Navy  (RE&S) 

Department  of  the  Navy 
Room  4E736,  Pentagon 

Washington,  D.C.  20350  1 

Commander,  Naval  Air  Systems  Command 
Attn:  Code  310C  (H.  Rosenvasser) 

Department  of  the  Navy 

Washington,  D.C.  20360  1 

Defense  Technical  Information  Center 

3uilding  5,  Cameron  Station 
Alexandria,  Virginia  22314  12 

Dr.  Fred  Saalfeld 

Chemistry  Division,  Code  6100 

Naval  Research  Laboratory 

Washington,  D.C.  20375  1 


No. 

Copies 

U«S.  Army  Research  Office 
Attn:  CRD-AA-IP 

P.0.  Box  1211 

Research  Triangle  Park,  N.C.  27709  1 

Naval  Ocean  Systems  Center 

A-ttn:  Mr.  Joe  McCartney 

San  Diego,  California  92152  1 

Naval  Weapons  Center 
Attn:  Dr.  A.  B.  Aoster, 

Chemistry  Division 

China  Lake,  California  93555  1 

Naval  Civil  Engineering  Laboratory 

Attn:  Dr.  R.  W.  Drisko 

Port  Hueneme,  Jalifornia  93401  1 

Department  of  Physics  &  Chemistry 

Naval  Postgraduate  School 

Monterey,  California  93940  1 

Dr.  A.  L.  S lafkosky 
Scientific  Advisor 
Commandant  of  the  Marine  Corps 
( Code  RD-i ) 

Washington,  D.C.  20330  I 

Office  of  Naval  Research 
Attn:  Dr.  Richard  S.  MjLller 

800  N.  Quincy  Street 

Arlington,  Virginia  22217  I 

Naval  Ship  Research  and  Development 
Center 

Attn:  Dr.  G.  Sosmajian,  Applied 
Chemistry  Division 

Annapolis,  Maryland  21401  1 

Naval  Ocean  Systems  Center 
Attn:  Dr.  3.  Yamamoto,  Marine 
Sciences  Division 

San  Diego,  California  91232  1 

Mr.  John  Doyle 

Materials  Branch 

Naval  Ship  Engineering  Center 

Philadelphia,  Pennsylvania  19112  1 


A 


SP4 72 -3/ A3 


472  :  GAN*.  716 : ddc 
78u472-608 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  GEN* 


No- 

Copies 


Dr-  Rudolph  J.  Marcus 
Office  of  Naval  Research 
Scientific  Liaison  Group 
American  Embassy 

APO  San  Francisco  96503  1 

Mr.  James  Kelley 
DTNSRDC  Code  2803 

Annapolis,  'iaryland  21402  1 


2 


SP4 72-3/B7 


472sGAKs7l6sddc 

78u472-608 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  356A 


No . 

Copies 


No. 

Copies 


Mr.  Robert  W.  Jones 
Advanced  Projects  Manager 
Hughes  Aircraft  Company 
Mail  Station  D  132 

Culver  City,  California  90230  1 

Dr.  C.  Giori 

IIT  Research  Institute 

10  West  35  Street 

Chicago,  Illinois  60616  1 

Dr.  M.  Litt 

Department  of  Macromolecular  Science 
Case  Western  Reserve  University 
Cleveland,  Ohio  44106  1 

DepartSfe^of  of  Science 

and  Meta^^!>g^^^Engineering 
Univer^i^of  Cinc^ftta^J. 

Cip^fnnati,  Ohio  45221*  1 

Dr.  Robert  E.  Cohen 
Chemical  Engineering  Department 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139  1 

Dr.  T.  P.  Conlon,  Jr.,  Code  3622 
Sandia  Laboratories 
Sandia  Corporation 

Albuquerque,  New  Mexico  1 

Dr.  Martin  Kaufmann,  Head 
Materials  Research  Branch,  Code  4542 
Naval  Weapons  Center 

China  Lake,  California  93555  1 

Professor  S.  Senturia 
Department  of  Electrical  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139  1 


Dr.  T.  J.  Reinhart,  Jr.,  Chief 
Composite  and  Fibrous  Materials  Branch 
Nonmetallic  Materials  Division 
Department  of  the  Air  Force 
Air  Force  Materials  Laboratory  (AFSC) 
Wright-Patterson  AFB,  Ohio  45433  1 

Dr.  J.  Lando 

Department  of  Macromolecular  Science 
Case  Western  Reserve  University 
Cleveland,  Ohio  44106  1 

Dr.  J.  White 

Chemical  and  Metallurgical  Engineering 
University  of  Tennessee 

Knoxville,  Tennessee  37916  1 

Dr.  J.  A.  Manson 
Materials  Research  Center 
Lehigh  University 

Bethlehem,  Pennsylvania  18015  1 

Dr.  R.  F.  Helmreich 
Contract  RD&E 
Dow  Chemical  Co. 

Midland,  Michigan  48640  1 

Dr.  R.  S.  Porter 
Department  of  Polymer  Science 
and  Engineering 
University  of  Massachusetts 
Amherst,  Massachusetts  01002  1 

Professor  Garth  Wilkes 
Department  of  Chemical  Engineering 
Virginia  Polytechnic  Institute  and 
State  University 

Blacksburg,  Virginia  24061  1 

Dr.  Kurt  Baum 
Fluorochem  Inc. 

6233  North  Irwindale  Avenue 

Azuza,  California  91702  1 


Professor  C.  S.  Paik  Sung 
Department  of  Materials  Sciences  and 
Engineering  Room  8-109 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139  1 


SP472-3/B5 


472:GAN:7l6:ddc 
78u472-6 08 


[ 

i 


r 


i 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  356A 


No. 

Copies 

Dr.  Stephen  H.  Carr 

Department  of  Materials  Science 

Northwestern  University 

Evanston,  Illinois  60201  1 

Dr.  M.  Broadhurst 

Bulk  Properties  Section 

National  Bureau  of  Standards 

U.S.  Department  of  Commerce 

Washington,  D.C.  20234  2 

Professor  G.  Whitesides 
Department  of  Chemistry 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139  1 

Professor  J.  Wang 
Department  of  Chemistry 
University  of  Utah 

Salt  Lake  City,  Utah  84112  1 

Dr.  V.  Stannett 

Department  of  Chemical  Engineering 
North  Carolina  State  University 
Raleigh,  North  Carolina  27607  1 

Dr.  D.  R.  Uhlmann 
Department  of  Metallurgy 
and  Material  Science 
Massachusetts  Institute 
of  Technology 

Cambridge,  Massachusetts  02139  1 

Naval  Surface  Weapons  Center 
Attn:  Dr.  J.  M.  Augl, 

Dr.  B.  Hartman 
White  Oak 

Silver  Spring,  Maryland  20910  1 

Dr.  G.  Goodman 

Globe  Union  Incorporated 

5757  North  Green  Bay  Avenue 

Milwaukee,  Wisconsin  53201  1 

Professor  Ratsuo  Ishida 
Department  of  Macromolecular  Science 
Case-Western  Reserve  University 
Cleveland,  Ohio  44106  1 


No . 

Copies 

Picatinny  Arsenal 

Attn:  A.  M.  Anzalone,  Building  3401 
SMUPA-FR-M-D 

Dover,  New  Jersey  07801  1 

•Dr.  J.  K.  Gillham 
Department  of  Chemistry 
Princeton  University 

Princeton,  New  Jersey  08540  1 

Douglas  Aircraft  Co. 

Attn:  Technical  Library 

Cl  290/36-84 
AUTO-Sutton 

3855  Lakewood  Boulevard 

Long  Beach,  California  90846  1 

Dr.  E.  Baer 

Department  of  Macromolecular 
Science 

Case  Western  Reserve  University 
Cleveland,  Ohio  44106  1 

Dr.  K.  D.  Pae 

Department  of  Mechanics  and 
Materials  Science 
Rutgers  University 

New  Brunswick,  New  Jersey  08903  1 

NASA-Lewis  Research  Center 

Attn:  Dr.  T.  T.  Serofini,  MS-49-1  1 

21000  Brookpark  Road 

Cleveland,  Ohio  44135 

Dr.  Charles  H.  Sherman 
Code  TD  121 

Naval  Underwater  Systems  Center 

New  London,  Connecticut  1 

Dr.  William  Risen 
Department  of  Chemistry 
Brown  University 

Providence,  Rhode  Island  02192  1 

Dr.  Alan  Gent 
Department  of  Physics 
University  of  Akron 

Akron,  Ohio  44304  1 


